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Resumen: 

Se investiga, mediante modelos matemáticos, los procesos de transporte de masa en una solución de 
muchos componentes. La solución se encuentra ubicada en una cámara limitada por dos barreras 
semipermeables, proximal y distal. A través de esas barreras se impone un flujo constante de sol- 
vente que arrastra los solutos hacia la barrera distal y produce contra flujos difusivos. Se define una 
región de extracción adyacente a la barrera distal y se establecen una desigualdad, relacionada con 
el comportamiento temporal del número de moles de los solutos en la región de extracción, tal que 
si se cumple asegura que la función de enriquecimiento relativo presente una fase transitoria por 
exceso respecto de su valor estacionario final. Teniendo en cuenta efectos de difusión cruzados, se 
construyen dos modelos matemáticos de los procesos de convección-difusión: un modelo tridimen- 
sional y un modelo global unidimensional. Se linealizan las ecuaciones de transporte y se resuelven 
en el marco de ambos modelos. Cuando los efectos cruzados son débiles, las soluciones de las 
ecuaciones de transporte de masa se simplifican mediante un procedimiento de perturbaciones 
regulares. A partir de las soluciones de las ecuaciones de ambos modelos matemáticos se investigan 
las condiciones en las que cabe esperar se verifique la mencionada desigualdad, que asegura que 
una extracción en un instante adecuadamente seleccionado del transitorio mejora el enriquecimiento 
relativo respecto de su valor en estado estacionario. 


Palabras clave: modelos matemáticos, procesos de transporte de masa, enriquecimiento de solu- 
ciones, teoría de perturbaciones, efectos de difusión cruzados, sistemas de ecuaciones a derivadas 
parciales, análisis modal. 


Abstract: 

Mathematical models are used to investigate mass transport processes in a solution of many compo- 
nents. The solution is in a chamber limited by two semipermeable barriers, proximal and distal. 
Through these barriers, a constant flow of solvent is imposed, which drags the solutes towards the 
distal barrier and produces diffusive counter flows. An extraction region adjacent to the distal barri- 
er is defined and an inequality is established, related to the temporal behavior of the number of 


moles of the solutes in the extraction region, such that if it is met, it ensures that the relative en- 
richment function presents a transient phase with an overshoot relative to its final stationary value. 
Taking into due account cross-diffusion effects, two mathematical models of convection-diffusion 
processes are constructed: a three-dimensional model and a one-dimensional global model. The 
transport equations are linearized and solved in the framework of both models. The solutions of the 
mass transport equations are simplified by a procedure of regular perturbations when the cross 
effects are weak. Based on the solutions of the equations of both mathematical models, the condi- 
tions which ensures that an extraction at an appropriately selected moment of the transient improves 
the relative enrichment with respect to its steady state value are investigated. 


Key words: mathematical models, mass transport processes, solution enrichment, perturbation 
theory, cross diffusion effects, systems of partial differential equations, modal analysis. 


(1) Introducción 


Consideremos un tubo cilíndrico de paredes impermeables y supongamos que en su inte- 
rior dos barreras semipermeables transversales al eje del tubo delimitan una cámara, que 
denominaremos cámara de enriquecimiento (Figura 1). 


extremo extremo 
proximal distal 


Figura 1 Esquema de una cámara de enriquecimiento por convección-difusión de longitud / y 
sección transversal circular de área S . En los extremos proximal (situado en x= 0 ) y distal 


(situado en x=1) respecto del flujo de solvente se ubican dos barreras semipermeables, 
atravesadas por un flujo constante de solvente desde el inicio del proceso. 


En dicha cámara se encuentra una solución con más de un soluto disuelto en un solvente. 
A través del sistema se hace atravesar un flujo de volumen constante del mencionado 
solvente, que penetra a través de la barrera semipermeable proximal y lo abandona a tra- 


vés de la barrera semipermeable distal. Asumimos que los solutos no pueden atravesar las 
mencionadas barreras. 

Supongamos que en el instante inicial la concentración de los diferentes solutos no varía 
de un punto a otro de la cámara de enriquecimiento: su distribución espacial es uniforme. 
Entonces el arrastre producido por el flujo tiende a aumentar las concentraciones de los 
solutos en las adyacencias de la barrera distal a expensas de una depleción de solutos en 
las adyacencias de la barrera proximal. 


Los gradientes de concentración de solutos que se producen en el interior de la cámara de 
enriquecimiento deben provocar a su vez contraflujos dispersivos que se oponen al flujo 
de convección, cada vez con mayor intensidad hasta que, asintóticamente, se establece un 
estado estacionario. 

En ese estado estacionario, para cada sustancia disuelta el flujo de arrastre por convección 
y el contraflujo dispersivo se anulan entre sí. 

En esas circunstancias cabe esperar que, en las proximidades de la barrera distal, una 
sustancia con menor coeficiente de dispersión se concentre relativamente más que otra 
con mayor coeficiente de dispersión. 

Como consecuencia, en estado estacionario, una extracción de una porción de solución 
adyacente a la barrera distal suministraría una solución enriquecida respecto de los solu- 
tos de menor coeficiente de dispersión. 

Sí en condiciones adecuadas esto efectivamente ocurre así, permitiría implementar un 
proceso de enriquecimiento de soluciones por efecto de los gradientes diferenciales de 
concentración en estado estacionario. 


Ahora bien, en igualdad de las demás condiciones, la velocidad con la que el campo de 
concentraciones de un soluto tiende a la distribución estacionaria correspondiente depen- 
de del coeficiente de dispersión de ese soluto y eventualmente de los efectos cruzados 
entre el soluto considerado y los demás solutos. 


Supongamos que fuera posible controlar los parámetros del proceso de modo tal que un 
soluto que se concentrara más en las proximidades de la barrera distal se aproximara más 
rápido a su estado estacionario. 


Entonces efectuando una extracción en la zona adyacente a la barrera distal en un instante 
adecuadamente elegido durante la fase transitoria del proceso se podría obtener una solu- 
ción aún más enriquecida en este componente que la solución obtenida en una extracción 
en estado estacionario (Figura 2). 
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Figura 2. Cociente de concentraciones de dos solutos (soluto 1 y soluto 2) en un punto situado 
en el interior de la cámara de enriquecimiento, lo bastante próximo a la barrera distal. Se ha 
representado una fase transitoria por exceso cuyo valor máximo se alcanza en el instante de 


tiempo í,, a partir del inicio del proceso. 


Esto constituiría una modificación no trivial del proceso original de enriquecimiento de 
soluciones, que podría denominarse proceso de enriquecimiento por efecto combinado de 
gradientes de concentración y control de los tiempos de relajación. 


La posibilidad de que un soluto que se concentra más en las proximidades de una barrera 
semipermeable pueda aproximarse más rápido a su estado estacionario se le ocurrió al 
presente autor en 1987, mientras revisaba las soluciones de un modelo matemático de un 
proceso de transporte de masa que había hecho en la década de los setenta. Este trabajo, 
de índole teórica, no fue publicado. Se llevó a cabo con la intención de comprender algu- 
nos de los resultados de unos experimentos sobre electroósmosis en células vegetales 
efectuados en 1970 por L. M. Sánchez y el autor. 


El desarrollo preliminar de la propuesta de modificar el proceso de enriquecimiento se 
hizo en un informe presentado por el autor a la Comisión Nacional de Energía Atómica de 
Uruguay [1] y en dos informes de investigación posteriores, más detallados y en coautoría 
con tres compañeros de trabajo en la Dirección Nacional de Tecnología Nuclear del Mi- 
nisterio de Industria y Energía de Uruguay [2] [3]. En la referencia [4] (páginas 84 y 85) 
se pueden encontrar una breve historia del desarrollo de los trabajos preliminares mencio- 
nados y una anécdota relacionada con una repercusión inesperada, en el campo nuclear, 
que tuvo la idea. 

Tanto en el informe original como en los informes más detallados que le siguieron, no se 
tuvieron en cuenta posibles efectos de difusión cruzados entre el soluto considerado y los 
demás solutos. 


El presente trabajo posee dos objetivos principales. 

Uno de ellos consiste, a partir de la introducción de una región de extracción adyacente a 
la barrera distal, en establecer una desigualdad, relacionada con el comportamiento tem- 
poral del número de moles de los solutos en esa región de extracción, tal que si se cumple 
asegura que la función de enriquecimiento relativo presente una fase transitoria por exce- 
so respecto de su valor estacionario final. 

El otro objetivo principal consiste en construir modelos matemáticos de los procesos de 
transporte de masa en la cámara de enriquecimiento para tener en cuenta los efectos cru- 
zados. Esto permitirá profundizar el análisis fisicomatemático de las condiciones bajo las 
cuales cabe esperar que el proceso modificado conduzca a un incremento en el enrique- 
cimiento respecto del que se obtendría trabajando en estado estacionario con la misma 
solución. 

En las referencias [2] y [3], además de un abordaje analítico, se emplearon diversas técni- 
cas de simulación digital para explorar las propiedades de las soluciones de los procesos 
de transporte de los solutos en ausencia de acoplamiento difusivo. En el presente artículo 
la investigación de las propiedades de las soluciones de los modelos matemáticos se abor- 
dará exclusivamente mediante herramientas analíticas, quedando el análisis numérico y 
las corridas de simulación digital para trabajos futuros. 


(2) Condición para que el enriquecimiento relativo presente una fase 
transitoria por exceso 


Para comenzar, es conveniente definir un parámetro de enriquecimiento adecuado para 
discutir el proceso en la cámara de convección-difusión. 

Con el propósito de conectar la discusión que se lleva a cabo en este trabajo con la inves- 
tigación del proceso de enriquecimiento tal como se llevó a cabo en los trabajos corres- 
pondientes a las referencias [2] y [3] y poder hacer referencia directa a los resultados 
obtenidos en los trabajos mencionados, comenzaremos considerando una solución terna- 
ria con dos solutos cuyas concentraciones (molares) uniformes iniciales son C,, y Cp, - 


Representemos por c, (r, F ) y C, (r, F ) las concentraciones (molares) de esos mismos 
solutos en un instante £ (medido a partir del inicio del proceso) en un punto de la cámara 
de vector de posición Y . Entonces definimos para el instante £ y el punto Y la función de 


enriquecimiento relativo (%,, del soluto 1 respecto del soluto 2 así: 


E A A, Q.1) 


aalt,7)= CANO 


Puesto que C; (0, E ) =C0 Y C, (o, F ) = Co, de la definición (2.1) resulta, en todo punto 


de la cámara de enriquecimiento: a,[0, F ) =1 Q.2) 


c,(t,7) c,(t,7) 


Siempre que —7—=X sea mayor que 53, O sea siempre que 0, (t,7)> 1, se 
c,(0,7) c,(0,7) 


puede decir que en ese punto y en ese instante la solución se encuentra enriquecida en el 
soluto 1 respecto del soluto 2, tomando como referencia la distribución uniforme inicial. 


A su vez, cuanto mayor sea 0%, lt E ), tanto mayor va a ser el enriquecimiento en el 


soluto 1. 
Cuando se alcanza un estado estacionario, la función de enriquecimiento en el punto F 
de la solución toma el valor correspondiente a las concentraciones estacionarias locales 


0 (7 ) Y Cao (7 ) de los solutos 1 y 2, respectivamente: 


a Oo EE 

Oro (7) ==. uo) Q.3) 
Cor Co 6 ) 

Para un flujo de volumen constante a través de la cámara, instalado a partir del instante 

inicial, Qt, (4 ed ) comienza siempre tomando el valor 1 y se aproxima asintóticamente a 


Arco (7). Cuando (%;) ., (7 )> 1, se puede decir que en estado estacionario la solución 
en el punto considerado se encuentra enriquecida en el soluto 1 respecto del estado de la 
solución al comienzo del proceso. Si (%,, ¿(E )<1 ocurre lo contrario: en estado esta- 


cionario se encuentra enriquecida en el soluto 2. 


De la investigación de las propiedades de las soluciones de un modelo matemático global 
unidimensional' del proceso de convección-difusión (desarrollado en las referencias [1], 
[2] y [3] empleando métodos tanto analíticos como numéricos) se desprende que si el 
punto considerado, de abscisa X (medida sobre el eje de la cámara) se encuentra lo sufi- 
cientemente próximo a la barrera distal, si el soluto 1 difunde con mayor dificultad que el 
soluto 2, y si se desprecian posibles efectos cruzados: 


(a) 7. (x)>1 
(b) Si además los números de Péclet de ambos solutos son mayores que 27, a, (e 5 x) 
evoluciona desde su valor inicial 1 hasta su valor estacionario final (%,» ax) presentan- 


do una fase transitoria por exceso respecto del valor final: para un instante intermedio 
tf =1,, la función onll , x) pasa por un máximo antes de aproximarse, decreciendo mo- 


nótonamente, a su valor estacionario final. Como consecuencia, para ese instante de tiem- 
po £,, y en igualdad de las demás condiciones, el enriquecimiento es máximo. 


| En un modelo global unidimensional las concentraciones de los solutos se consideran dependien- 
tes de la posición a lo largo de un eje que conecta la barrera proximal con la barrera distal. Pueden 
interpretarse como promedios sobre las secciones transversales a este eje. 


En última instancia, para evaluar las posibilidades del proceso hay que diseñar, construir 
y Operar una unidad piloto, en la cual se deberá extraer un cierto volumen de una solución 
de muchos componentes en un instante y en un lugar adecuado. 

Representemos mediante R a la región espacial, fija al laboratorio y ocupada por la cá- 
mara de enriquecimiento de esa unidad piloto. Consideremos ahora una región de solu- 


ción R, fija respecto del laboratorio (contenida en R ) cuyo volumen v(R,) se extrae 
en el instante f. 


El número de moles », del soluto K presentes en esa región en el momento en el cual se 


efectúa la extracción viene dado por la integral de volumen extendida a la región R, : 


n, (t,R,)= fc, (t,7). av (2.4) 
Ro 
El número de moles en esa misma región al inicio del proceso viene dado por: 
n.(0;R,)= |c,(0,F)- dV =cy, «V(R,) (2.5) 
R 


e 


Para una solución de M solutos en un solvente, introduzcamos ahora los cocientes de 
concentración para cada soluto, referidos a la concentración uniforme inicial: 


0 (0,7) 17) (i=1,2,...,n) 2.6) 
Co¡ 
De (4), (5) y (6) se desprende: nltsR,) 5 J o lt F)-dV  Q7) 


n.(0;R,) V(R,) 


Entonces es posible definir la siguiente función de enriquecimiento relativo de la solución 


Ro, 


en el soluto K respecto del soluto [ en la región de extracción Ro: 


A o) fo. 1,F)-dV E 


al fol 1,F)-dV 


También en este caso se verifica (X,, (0; R.) = 1 mientras que: 


ae 
O. (R, )= lim 0, ( f; R,) av 
oo 


0%, (e ¿Ro % 


(2.9) 


La factibilidad de producir circunstancias en las cuales cabe esperar que 0, (t s R,) pre- 


sente una fase transitoria por exceso, constituye la base de la modificación propuesta al 
proceso de enriquecimiento en una cámara de convección-difusión. 


Para analizar esas circunstancias, comencemos por descomponer el número de moles de 
un soluto presentes en cada instante en R, de esta forma: 


nleR )=n, (AR) +Sn (ER) (¡=12,...,n) (2.10) 


Dadas las condiciones iniciales y la naturaleza del proceso de arrastre por convección, 
cabe esperar que n, 6 ¿R,) aumente progresivamente desde su valor inicial (2.5) hasta su 


valor estacionario final A, ¿(R.). Entonces Ó n; 6 R,)debe ser una función negativa y 


0% ¡R,) 
Oro (R.) 


res mayores que 1, entonces (%;, (e > R,) presenta una fase transitoria por exceso respecto 


monótona creciente con límite O para 1 >-+00, 


Si desde un instante de tiempo en adelante el cociente tiende a 1 desde valo- 


de su valor estacionario final, alcanzando un valor máximo para algún instante f =1,,. 
As lr > R.) 


Ol (R, ) 


. Pero como Sn, (1 ¿R,) y Sn, (t ¿R.)son funcio- 


Aplicando la descomposición (10) a los solutos K y Íl resulta que 


e Sn, (t,R,) 3 Sn, (t,R,) 
My. (R,) n,..(R,) 


nes siempre negativas, la última desigualdad equivale a la siguiente: 


Sn, (tR,) Sn,(t;R,) 
< 
Mos (R,) m,. (R,) 


Si desde un instante de tiempo en adelante se verifica (2.11), el coeficiente de enriqueci- 


Q.11) 


miento Q salt ¿R,) presentará una fase transitoria por exceso. 


Para poder aplicar esta condición en el marco de una investigación teórica, es necesario 
calcular (2.10) a partir de un modelo matemático del proceso en la cámara de enriqueci- 
miento. 


(3) Construcción de modelos matemáticos del proceso de enriquecimien- 
to que tienen en cuenta los efectos de difusión cruzados 
(3.1) Modelo espacial tridimensional 


Consideremos una solución con 1 +1 componentes, N solutos y el solvente. Las concen- 
traciones molares de los solutos, en un instante f (medido a partir del inicio del proceso) 


en un punto de la cámara de vector de posición Y, vienen representadas por C, (e mA ) para 


¡=1,2,...,n. Para cada soluto, la conservación de la masa, en ausencia de reacciones 


químicas, viene dada por [5] [6] ná =V ey, ¡=12,....n (G.1) 
Í 


En esta ecuación J : (e cd ) representa el campo vectorial de densidad de flujo del soluto 


respecto del laboratorio, debida a la combinación de arrastre por convección y difusión y 
V e es el operador divergencia [5] [6]. 


Para describir la difusión en una solución de muchos componentes se pueden utilizar tres 
sistemas de referencia, según como se obtenga la velocidad de convección: promediando 
sobre los moles, promediando sobre la masa o promediando sobre el volumen [7] [8]. Los 
valores que se asignan a los coeficientes de difusión dependen del sistema de referencia 
adoptado. Se dispone de ecuaciones que permiten el cálculo de los coeficientes de difu- 
sión respecto de un sistema de referencia en función de los coeficientes de difusión res- 
pecto a otro sistema [7] [9]. 

En general las mediciones de los coeficientes de difusión se llevan a cabo referidas a un 
sistema fijo al volumen [10] [11]. Por esta razón y porque la matriz de coeficientes de 
difusión (que se introduce y discute más adelante) cuando se la refiere al volumen es más 
fácil de analizar, en el presente trabajo se emplea el sistema fijo al volumen para describir 
los flujos de difusión. 

La velocidad con la que se produce el arrastre por convección de la solución referida al 


volumen se representa por ve ,F ) mientras que la densidad de flujo difusivo del soluto 


La oa Vr 
respecto del movimiento por convección viene dada por J, (e z 7) 3 


La densidad de flujo respecto del laboratorio debido a la combinación de arrastre por 
convección y difusión del soluto resulta entonces: 


a po E 

J (e,F)=c,(c,F)-v(t,7)+J7,(t,F) (3.2) 
Asumiendo que la cámara de enriquecimiento se comporta como un sistema isotérmico e 
isobárico en el cual la difusión se debe a los gradientes de concentración, las densidades 
de flujo de los solutos se pueden relacionar con sus concentraciones mediante una genera- 


lización de la primera ley de Fick, aplicable tanto a soluciones líquidas como a mezclas 
gaseosas [10] [11] [12] [13] [14] [15]: 


Til t,F)= 0 (c,(1,5),....c, (t,7))- Ve, (t,F) (3.3) 


Los coeficientes D;, se suelen denominar coeficientes de difusión principales, mientras 
es Vv . . . dia . ., 
que los coeficientes D; con 1% /se denominan coeficientes de difusión cruzados. 


En las paredes de la cámara, incluyendo las barreras semipermeables proximal y distal, 
los flujos de los solutos deben ser nulos, debido a que esas paredes (se supone) son im- 


permeables para los solutos: J, (r, 7)o n=0 ( =D n) (3.4) 


En esta condición en la frontera M es un vector unitario (versor) perpendicular al plano 
tangente al punto de la frontera considerado?. Las concentraciones iniciales son constantes 
en todo el interior de la cámara: 


c/(0,7)=c,, (i=1,2,...,n) (3.5) 


Con el fin de poder efectuar un abordaje analítico, en lo que sigue utilizaremos los valores 
de los coeficientes de difusión para las concentraciones iniciales de los solutos: 


Diaria AISD. (Cuscs.) (i=1,2,....n)  G.) 


a Ed A 5 v 
De esta forma los elementos de la matriz de coeficientes de difusión D = [D; ] son 
nxn 


constantes. Esta simplificación ha sido utilizada previamente con buenos resultados [10] 
[11]. 
Asumimos que el fluido es incompresible: Ve vr E ) =0 (3.7) 


El solvente no puede atravesar la pared lateral de la cámara, debido a lo cual, si M es un 
versor normal a la cara interna de la pared en el punto considerado, la componente normal 


de la velocidad en ese punto debe anularse: ve, 7)o n=0 (3.8) 
Si la solución moja la pared, la velocidad misma debe anularse allí. 
Para hallar el campo de velocidades ve ¿E ) se debe establecer una conexión con las 


ecuaciones de la mecánica de fluidos. Para establece esa conexión conviene describir la 
difusión mediante un sistema de referencia fijo al centro de masas en lugar de un sistema 
de referencia fijo al volumen como el utilizado en este trabajo. Pero si se va a resolver el 
problema por simulación digital, empleando un código de cálculo, se puede efectuar una 
transformación de sistema de referencia como se mencionó previamente. 


A partir de las relaciones (3.1), (3.2), G.3), (3.6) y 3.7) se obtienen las ecuaciones que 
describen la evolución de las concentraciones de soluto en la cámara de enriquecimiento: 


7 (1,7) +v(1,7)e Ve,(+,7)= YD «We, (e,r) 69) 
jal 
Si se supone que ve E ) viene dado de antemano, las ecuaciones de evolución (3.9) junto 


con las condiciones en la frontera (3.4) y las condiciones iniciales (3.5), dan origen a un 
modelo matemático en tres dimensiones espaciales, que se puede utilizar para describir 
los procesos de transporte de solutos en una cámara de enriquecimiento, teniendo en 
cuenta los efectos de difusión cruzados. 

En principio este modelo puede aplicarse a cámaras de enriquecimiento arbitrarias (no 
necesariamente cilindros) y a flujos de solución tanto laminares como turbulentos. 

En lo que sigue el campo de velocidades de transporte de volumen se asume es estaciona- 


rio, correspondiendo a un flujo laminar: ve E ) = v(7 ) (3.10) 


2 En las referencias [16] y [17] se discute y justifica este tipo de condiciones de borde para modelos 
matemáticos de procesos de transporte de masa. 


O a ZA pvp Ss 
En ese caso: ay ele r)+v(5)o Ve (1,F)=2.D, -V c,(t,F) (3.10 
j=1 
No obstante, parece que aun en el caso de una cámara de enriquecimiento con forma de 
cilindro y asumiendo la máxima simetría posible para este caso, el análisis de la fase 


transitoria de la función de enriquecimiento 0, (e ¿R,) a partir de las soluciones analíti- 


cas para las ecuaciones (3.11) solo puede llevarse a cabo en forma parcial. Para obtener 
información más detallada de lo que acontece es necesario-efectuar un abordaje mediante 
técnicas de simulación digital. 


Para poder llevar más lejos el abordaje analítico, como se pretende en el presente artículo, 
es necesario simplificar el modelo matemático del proceso de enriquecimiento. 


Una simplificación posible consiste en promediar? los campos (de velocidad v(F ) y de 


concentraciones de soluto C, 6 , 7) sobre las secciones transversales al eje de la cámara 


de enriquecimiento, que de ahora en adelante se supone tiene la forma de un cilindro. Con 
la operación de promediado aplicada al modelo tridimensional se puede obtener un mode- 
lo global unidimensional mucho más simple para analizar, pero que permite tener en 
cuenta los efectos de difusión cruzados. 


(3.2) Modelo global unidimensional a partir del modelo tridimensional 


Para una cámara de enriquecimiento con forma de cilindro de radio a y longitud l, en 
condiciones de flujo laminar, la magnitud del campo de velocidades puede representarse 


2 
mediante la fórmula bien conocida [6]: v(r) = Deal = (3.12) 


En (3.12) la distancia del punto considerado al eje del cilindro se representa por / mien- 
tras que y representa la velocidad promedio del flujo de volumen a través de una sección 
transversal al eje de la cámara. El campo de velocidades locales es independiente de la 


posición X (comprendida entre x=0y x=1) de la sección transversal respecto del eje 
del cilindro. 


Si e. es un versor a lo largo del eje y e, es el versor radial correspondiente al punto con- 


siderado, las densidades de flujo de soluto se reducen a: 


Ji (e, r,x) = (elena En: le napa Eo elena) 
5 Ox Ga Or 


(3.13) 
Combinando la conservación de la masa (3.1) y con las fórmulas para las densidades de 
flujo (3.13), se deduce las ecuaciones para la evolución de las concentraciones: 


3 En la referencia [5] se puede hallar un desarrollo exhaustivo de los métodos de promediado apli- 
cados al estudio de los procesos de transporte de masa. 


otr) 01) e ltsria)= YD) (Lrzot Ds (r, 0) (3.14) 


La condición en la frontera (3.4) se reduce a J P (e , 0, x)e e = 0 sobre la cara curva de 


la cámara. 


Teniendo en cuenta (3.13), esto equivale a: 20 ¿(e a, x)= 0 (3.15) 


Como se puede suponer que la matriz de E de difusión es invertible [], de 


(3.15) resulta para todo 1 > 0 y todo x: Ze, (t,a, x) =0 li =1, 2,...,N) (3.16) 
r 


Dado un campo X 6 Y, x) (que puede representar concentraciones, velocidad, densidad 


de flujo) construimos su promedio sobre una sección transversal al eje de la cámara: 


X(t,x)= [Xle,r,x)-270r-dr (3.17) 


a? 
Denominemos ct, x) ¿E (e ,x) y V alas concentraciones promedio, densidades de 


flujo promedio y velocidad de flujo de volumen promedio. 

Promediando los flujos a partir de las ecuaciones (3.13), admitiendo que por simetría se 
anula el promedio de la componente radial de la densidad de flujo y aplicando la teoría de 
Taylor y Aris de la dispersión por difusión molecular en un flujo laminar de solución 


[16][18][19] se obtiene para el promedio q Ke, x) de la componente longitudinal del 


flujo: 7 (t,x)=G (t,x)-5-D", - z a (t, x) Y Dr. 2, (t, x) (1=12,....1) 
x 


ji, j=l 

(3.18) 

El coeficiente principal efectivo de difusión Dei , £eneralmente denominado coeficiente 
de dispersión, se expresa en términos de parámetros en principio conocidos previamente: 
a? y? 


D", =D] + pi (i=1,2,...,n) (3.19) 
Promediando la ecuación de conservación de la masa (2.1) se obtiene: 
O _ O 5 
€; =-/4, ¡=1,2,...,n (3.20) 
Ot Ox ) 


A partir de (3.18) y (3.20) se deduce finalmente la ecuación que describe la evolución de 
los campos de concentraciones a de los pe 


ón 0 j- 
a ex) +72 E(t,x) =D", - - >El da LE (ex) =1,2,....1)  G.21) 


Las condiciones en la frontera se reducen ahora a restricciones impuestas sobre las barre- 
ras proximal y distal, para todo 120:  J, (1,0) =J, (2,1) =0  (í=1,2,....n) (3.22) 


Las condiciones iniciales (3.5) se aplican a las concentraciones promediadas: 
G.(0,x)=c,, (i=1,2,....n) (3.23) 
Si se supone que Y viene dada de antemano, las ecuaciones de evolución (3.21) junto con 
las condiciones en la frontera (3.22) y las condiciones iniciales (3.23) constituyen un 
modelo matemático simplificado a una dimensión espacial (modelo global unidimensio- 
nal) que se puede utilizar para describir en forma resumida los procesos de transporte de 


solutos en la cámara de enriquecimiento, teniendo en cuenta los efectos de difusión cru- 
zados. 


Si se asume que los efectos de difusión cruzada son despreciables, cosa que puede acon- 
tecer en una solución muy diluida, el flujo promediado sobre una sección transversal al 
eje de la cámara se reduce a: 


pe 2 a O _ z 
J (e, x)= 3, (t,x)-v -D?, -— E (t, x) M2) (3.24) 
> OA 
La ecuación de evolución de la concentración promediada sobre la misma sección se 
reduce a: 


2 

E E A A (3.25) 
Ot Ox "Ox 

Las ecuaciones desacopladas (3.24) y (3.25) constituyen la base del modelo matemático 

estudiado en las referencias [1], [2] y [3]. 

En esos trabajos las ecuaciones (3.24) y (3.25) se introdujeron postulándolas ab initio, sin 

deducirlas a partir de un modelo tridimensional del transporte. Por este motivo el coefi- 


ciente principal efectivo de difusión D, ,, (coeficiente de dispersión) se introdujo direc- 


tamente como un parámetro fenomenológico sin relacionarlo con el coeficiente de difu- 
sión molecular, ni con el radio de la cámara, ni con la velocidad promedio de la solución. 


(4) La matriz de coeficientes de difusión. Problema de valores y vectores 
propios, modos normales de difusión binaria y soluciones de los mo- 
delos matemáticos del proceso de enriquecimiento. 


Para una solución de +1 componentes la matriz de coeficientes de difusión se puede 
representar así: D= [p, | (4.1) 
Y dnxn 


Esta matriz en general no es simétrica. 

Mientras que se dispone de mucha información experimental sobre los coeficientes de 
difusión en soluciones binarias, la medición de coeficientes de difusión en soluciones de 
más de dos componentes se ha limitado a soluciones ternarias y cuaternarias, y aún en 
estos casos la información originada en experimentos es escasa [11]. 

Por otra parte, el desarrollo de la teoría no ha sido todo lo detallado que cabría esperar, 
teniendo en cuenta la importancia de los efectos de difusión cruzados en procesos indus- 
triales y en procesos ambientales [20]. 


Consideremos el problema de valores y vectores propios para la matriz de coeficientes de 
difusión: D. U, = D, -U, (4.2) 


Cuando se utiliza un marco de referencia fijo al volumen y se asume que tanto la matriz 


- 


L =|1, | de coeficientes fenomenológicos de Onsager [7] [21] como la matriz 


- 0 


H=|—u, de las derivadas parciales de los potenciales químicos respecto de los 


On, 


J nxn 


números de moles [13] [14] son ambas definidas positivas*, los valores propios D, resul- 


tan reales y positivos [15]. 
En lo que sigue supondremos que esos valores propios son distintos dos a dos, de modo 


que el conjunto de vectores propios columna (ú,,ú,..... 8, ) es linealmente independien- 


¿ ; ; PT : 
te []. Además, los vectores propios se toman normalizados: si U n es el vector (fila) tras- 


=T - NE 
puesto 4, 4, =1. 


Introduzcamos ahora un vector columna para cada uno de los solutos: 


e, =col(0,....0,1,0,....0)=col(5,), , (i=12,....n) 


IS 
El símbolo O; representa la delta de Kronecker (0, =1, O; =0siiAH ]J). 
La proyección de cada uno de los vectores propios sobre cada uno de los vectores corres- 


pondientes a los n solutos origina una matriz invertible U' de nxn: 


U = [7 4 U, l.< lo, les (4.3) 
Consideremos ahora la posibilidad de construir una solución tentativa (ansatz) para los 
campos de concentraciones de los modelos tridimensional y global unidimensional a 
partir de ciertos campos auxiliares asociados a los valores propios y a los vectores propios 
de la matriz de coeficientes de difusión. 

Hay varios antecedentes del empleo de este enfoque en investigaciones de procesos de 
difusión molecular en soluciones de muchos componentes, en estado sólido y en líquidos 
[10] [11] [12] [14], pero al parecer este procedimiento no se ha desarrollado para el tipo 
de procesos de convección-difusión como los que cabe esperar ocurran en la cámara de 


La matriz simétrica de coeficientes fenomenológicos L (que relacionan los flujos termodinámicos 
con las fuerzas termodinámicas) es siempre definida positiva [7]. No obstante, parece que para 
ciertos estados termodinámicamente estables de algunas soluciones muy concentradas la matriz 


simétrica A puede no ser definida positiva. Entonces algunos de los coeficientes de difusión prin- 
cipales D, pueden ser negativos [11]. Estos casos serán excluidos en este trabajo. Cuando se em- 
plea un sistema de referencia fijo al volumen, D=L-S siendo S =P” -H-P con P inverti- 


ble. Si H es definida positiva, los valores propios de D son reales y positivos [14]. 


enriquecimiento. El análisis mediante campos auxiliares para la cámara de enriquecimien- 
to se desarrolla en esta parte (4). Se construyen soluciones analíticas para los dos modelos 
presentados en la parte (3): la solución del modelo tridimensional se construye en la sec- 
ción 4.1 y la solución del modelo unidimensional se construye en la sección 4.2. 


... 


(4.1) Solución para el modelo tridimensional y condición para que el coeficiente de 
enriquecimiento relativo presenta una fase transitoria por exceso respecto de su 
valor estacionario final. 


Consideremos entonces una posible solución tentativa (ansatz) para los campos de con- 
centraciones del modelo tridimensional a partir de ciertos campos auxiliares Pp, 6 ,F) 


asociados a los valores propios D A los vectores propios U A de la matriz de coeficien- 


tes de difusión: 
c (t,7)=é? - Y 9,(1,7)-5, a Ye (17) 0, (== 221) (4.4) 
ul p=l 


Sustituyendo el ansatz (4.4) en las ecuaciones (3.11) se obtienen las siguientes M ecua- 


ciones desacopladas para los campos auxiliares bp, (e A ) : 
9) NT ER El 
oyo e Fr) Va, lt,7)=D, -V*9, (1,5) (u =1,2,...,n) (4.5) 


Cada campo 0, 6 , 7) se comporta como lo haría la concentración de un soluto en una 
solución binaria, con un coeficiente de difusión D A sometido a un proceso de convec- 


ción caracterizado por el campo de velocidades v(7). Tomados colectivamente, estos 
campos auxiliares se pueden denominar modos normales de difusión. 


De la forma (4.4) asumida para los campos de concentraciones de solutos se desprende 
que todos esos campos se expresan como combinaciones lineales de los mismos campos 


auxiliares. La diferencia entre un soluto K y otro soluto / aparece en la diferencia que 


A T 
=U ,, ; y las proyecciones €, :u 


ÑO : >T > 
pueda existir entre las proyecciones €, :U a 


y = Ue 


Teniendo en cuenta las ecuaciones (3.2) y (3.3), la densidad de flujo del componente ¡ se 
puede expresar como combinación lineal de campos vectoriales de densidades normales 


de flujo Y, (t, F): T(0,7)= E (0,7), =0  (i=12...1) (4.6) 
u=1 


En (4.6), por definición: 
Y, (t,)=0,(c,F)-v(+)-D, -V4, (er) (u=12,....n) (4.7) 


Las restricciones sobre las densidades de flujo (ecuaciones (3.4)) que se imponen en la 
frontera OR de la región R equivalen a las siguientes condiciones impuestas sobre las 


densidades normales de flujo, para todo punto F perteneciente a OR: 


Y, (t,F)eñ=0 (u=1,2,....n) (4.8) 


a a O E OS 
Teniendo en cuenta (4.7), si V, (7) = v(7)e n y on b, (r, F) = VÓ, (r, F)o n resulta 
n 
Y se 9) E 
que (4.8) se puede poner así p,(t,F)- v, (7) = D, ] dalt, 7) (en OR) (4.9) 
n 

La solución de la ecuación (4.5) se puede expresar como suma de una solución estaciona- 
da (7) y un transitorio ó0, (e, r) que tiende a cero para l[ >+00. 
La solución estacionaria verifica la ecuación: 

v(r)o VE,..[(F)=D, -V*9,..(5) (4.10) 


El transitorio verifica la ecuación: 
0 O E ze 
Si 59, (t,F)+v(5)0V 89, (t,F)=D, -V?SH,(t,F) (4.11) 


Separando variables en esta última ecuación, buscando soluciones de la forma 


T 6 ) Y, (7 ) se obtienen estas dos nuevas ecuaciones: 


(0) =—0% -T, (1) (4.11) 
o, w,(r)+v(r)o Vy,(F)=D, Vw, (7) (en Ry (4.12) 
A las soluciones de (4.4) se les imponen la restricción v, (5 ). yw, (5 )= D,: - y, (5 ) 


sobre la frontera OR. Se obtiene así un problema de Sturm-Liouville que solo posee 


-z . 2 
solución para un conjunto numerable de valores de (Y, que se pueden ordenar como una 


sucesión creciente o, ¡< o, 9 <n. < o, e <... con límite infinito [22] [23]. 


Poniendo (0, , =7— resulta finalmente la siguiente solución como suma de una serie 


4,p 


de funciones: p,(t,7)= 0, (5) +00, (1,7)= Eo A (7) (4.13) 


Teniendo en cuenta (4.4), los campos de concentraciones de los solutos vienen dados 
por: 


cle,7)= 0, MAD: y, (a (i=1 2,1) (4.14) 


u=1 p=1 
En esta expresión el conjunto de tiempos de relajación T, p Para cada modo normal 


forma una sucesión monótona decreciente cuyo límite es 0: 
yd IE o E O E el (4.15) 
Asintóticamente: lim C, (t,7) r Y, (7 Es da (i = 1,2,...,n) (4.16) 


Pero como la cantidad de cada soluto en la cámara de enriquecimiento permanece inva- 
riante durante el proceso, se verifican las igualdades: 


[c,(t,7). dV=V(R)-c, (i=1,2....1) (4.17) 
R 
En particular: lima [a 1,F) )- dv = ¡E )- dv (í=1,2,...,n) 
R 


Entonces, de estas didas ecuaciones y de las ecuaciones (4.16) se desprende que las 


integrales de volumen de los valores asintóticos de los campos auxiliares J Pp, ¿(5)- dV 
R 
verifican el sistema de las ecuaciones: 


2 J9..0 Vu, =V(R)-c,, (6=1,2,.011) (4.18) 
u=1 
Como la matriz de coeficientes U ni es invertible, de (4.17) se obtienen las integrales 


Y 
> 


o, (5): dV (lu = 1,2.) en función de los valores iniciales de las concentraciones 
de los solutos en la cámara de enriquecimiento. 


A su vez, de (4.17) y (4.14), teniendo en cuenta (4.18), se deducen las siguientes restric- 
ciones sobre las funciones Y, , r 


[wa ,(7) )dV=0  (u=1,2..n) (p=12..) (4.19) 


En las fórmulas (4.14) el tiempo de relajación dominante Toa para cada modo normal 


determina el comportamiento asintótico de las concentraciones. Entonces, teniendo en 
cuenta (4.16), se tiene la siguiente aproximación asintótica para las concentraciones: 


t 
n == 


clt,7)=c, (FAY 3e wr, (5) pa, (i=1,2,....n) (4.20) 


u=l 


A partir de las ecuaciones (4.20) y (2.10) se pueden estimar las evoluciones temporales de 
los números de moles de los solutos en la región de extracción: 


n(eR,)=0n, (R)+Sn (ER) (i=1,2,....n) 
Los números de moles en estado estacionario se pueden calcular integrando los campos 
estacionarios sobre la región de extracción: Mn; ¿(R,) = J C; ¿(E ) dV (4.21) 


R 


e 


El comportamiento asintótico de los transitorios se puede estimar mediante las fórmulas: 


Sn(t;R,) =D” -Ju,lo» Ven, (i=1,2,....n) (4.22) 


H= 2, 


A partir de las fórmulas (4.21) y (422) se puede reescribir así la desigualdad (2.11): 


e Iv, tr) Vu, y e Ta Sy, .(5)-dV MU 
2, u=l R 
< E (4.23) 
A; o) (R, ) A; o (R, ) 


Si desde un instante de tiempo en adelante se verifica la desigualdad (4.23), el coeficiente 
de enriquecimiento relativo cl a R,) presenta una fase transitoria por exceso respecto 


Meco (R,) 


de su valor estacionario A. (R,)= —%— , como se mostró al final de la sección 
nel ) 


(2): aún en presencia de efectos de difusión cruzados el enriquecimiento relativo pasaría 
por un máximo durante el transitorio entre el estado inicial y el estado estacionario final. 


(4.2) Solución para el modelo global unidimensional 


Consideremos nuevamente la posibilidad de construir una solución tentativa (ansatz) para 
los campos de concentraciones, ahora para el modelo global unidimensional, a partir de 


campos auxiliares O, 6 , x) asociados a los valores propios D, y a los vectores propios 


17] > de la matriz de coeficientes de difusión: 


e [Entoa, Foca), (i=1,2,...,n) (424) 
u=1 u=1 


Sustituyendo el ansatz (4.21) en las ecuaciones (3.18) se obtienen las siguientes NM ecua- 


ciones desacopladas para los campos auxiliares Y, 6 ; x) : 


ai e Wei 485) 
Ot a 7 pee 


Teniendo en cuenta las ecuaciones (3.18), la densidad escalar de flujo J P 6 s x) del com- 
ponente i se puede expresar como combinación lineal de campos escalares de densidades 


o) 
=P lt,x): 


normales de flujo Y) = p,(t,x)-v =D: 0 
Xx 


)= t,x)-u,,=0  (i=12,...,n) (4.26) 


Las condiciones de borde en las barreras proximal (Xx = 0) y distal (x= L) se replantean 


como condiciones de frontera para los flujos auxiliares: a 6 0) = Ma 6 sl ) =0 (4.27) 


Introduzcamos las variables sin dimensiones, espacial £ y temporal 7 , mediante las 


l 
relaciones x=1-¿ y t= * *T, junto con el número de Péclet correspondiente a 
Vv 


cada modo normal de difusión: Pe = LA [7 A A n) (4.28) 

4 D, 
Entonces las ecuaciones (4.22) y las condiciones en la frontera se pueden reescribir res- 
pectivamente así: 


0 0 a 
nr E) ae E)= Pe, 0% 


d 
Pe, Ñ P, (7,0) + dé P, (7,0) Pe; Sl P, (7,1)= de 


La solución de la ecuación (4.29) se puede expresar como suma de una solución estacio- 
naria O... (x) y un transitorio Ó O, 6 , x) que tiende a cero para [ —>+00. 


(7,E) (u=12,....n) (4.29) 


Cole) (430) 


Los campos de concentraciones de los solutos, de acuerdo con (4.21) se pueden expresar 
así: 


E lt,x) +20, x)-u,, (4310) C,ola)=D O... (1), (4.32) 
.u=1 
De la ley de conservación de la masa a =-=J ¿ Se deprende que en estado esta- 
Í Xx 


cionario el flujo J, debe ser constante, independiente de la posición sobre el eje de la 


cámara. Como se anula en los extremos, debe ser nulo en todas partes. Entonces las solu- 


ciones estacionarias (,, (E) para los campos auxiliares deben verificar las ecuaciones: 


Pe, -9,.(E)= E E (E) (1 =1,2,...,n) (4.33) 
Las soluciones estacionarias correspondientes a las ecuaciones (4.33) son: 
0... ¿)=0,,..(0)-explPe, - ¿] (4.34) 


Para una cámara cilíndrica de longitud / las soluciones estacionarias C; (3) (promedia- 

das sobre las secciones transversales al eje) se relacionas con las concentraciones iniciales 
l 

a través de las fórmulas J C, (x): dx=1- Coj Sustituyendo (4.32) en esta última fór- 
0 

mula, y teniendo en cuenta (4.34) se obtiene el siguiente sistema de ecuaciones lineales 

para determinar los valores de Y, 0) lu =1, np 


yn - lu, Ou 0)=1-3,,  (i=12,...,n) (4.35) 
u=l Pe, 
Por su parte el transitorio Ó Ñe (r E ) verifica la ecuación: 
2) cs 0 
—Ó 20 «00, Lt 4.36 
Or P, la 5)+ OE p, (e ele Pe, Ox? P, (r $) ( ) 


Separando variables y buscando soluciones de la forma T, (c ) X de (€ ) se obtienen dos 


d 
relaciones. Para la función del tiempo resulta Ena T, (7) = -0, : de (z ) (4.37) 


e 2 
Su solución, a menos de una constante, es T, (c ) = exp|- 0,7 | (4.38) 


Para la función de la coordenada => se obtiene la ecuación: 


2.0%. 
2 E X= 5210) (O<É<I) (4.39) 


Las condiciones en la frontera (4.30) se reducen a las siguientes restricciones impuestas 
d d 

sobre X (€): Pe, X(0)=¿¿X 00) Pe, X= ¿XD (4.40) 

Si las soluciones de (4.39) deben verificar (4.40), los únicos valores admisibles para o, 


2 2 2 
_|Pe,+4-n"-z (4.41) 
4 Pe, 
2 


Pe, 0: X (£)+ Pe 


son los que forman la sucesión creciente: (o ) 
4,1 )n=1,2,... 


Teniendo esto último en cuenta, se obtiene la sucesión de soluciones: 


P . 
Xion (£)= Sl = j S sen|n -a ¿+ 0...) 186, ,]= = e (4.42) 
ly 


Combinando (4.38) con (4.42), teniendo en cuenta (4.41) y agregando una constante mul- 


, 2 cada término, el transitorio Sp,(r La ) en el modo normal se puede 


tiplicativa B, 
expresar como suma de una serie de funciones: 


Pe? +4n? 7? 
4 Pe, 


5p,(r,£)=Y B,, al 4-Pe,, E cena -senln-1-£+0,,] (4.43) 
n=1 


Las constantes B jon Se calculan a partir de la serie correspondiente al instante 7 =0: 
Pen: 


59, (0,)=0,(0,€)-0,.(£)=Y B,,e > «senno-£+0,,] 


n=1 


Los valores iniciales de los campos auxiliares Y, (0, 5 ) se determinan a partir del siste- 


ma de ecuaciones lineales que se obtiene de (4.24) poniendo £=0 e introduciendo la 


Xx n 
variable sin dimensiones ¿= y a (0,€) = Yo, (0, E)- Uni = Cos (4.44) 
u=1 
Se deducen finalmente las siguientes expresiones para las constantes: 
1 Pep:ó 
= Xx 
B,,=2:[5p,(0,É)-e ? - E “7 (>) + Ona -dé (4.45) 
0 


Entonces, de (4.24), (4.32) y (4.43), teniendo en cuenta las relaciones x=/-É¿ y 


l Ñ Ñ ; 
=(? . 7 , resulta que cada uno de los campos de concentraciones promediadas viene 
v 


dado por la siguiente fórmula: 
O ES 
ld o Ne (.). 20) senju() +0, 2u,, 446) 
u=l n=1 


Los tiempos de relajación de los modos normales vienen dados, para cada modo, por las 
sucesiones monótonas decrecientes a medida que M aumenta: 


E NO E E E sl 
PA A e A 


Las fórmulas (4.46) corresponden, para el modelo global unidimensional, a las fórmulas 
(4.14) del modelo tridimensional. Pero ahora se disponen de fórmulas analíticas para 


calcular los tiempos de relajación en función del índice H y del número de Péclet Pe ” 


correspondiente a cada modo normal. Examinando el comportamiento de T, , en función 


de Pe > Para A fijo, resulta que si incrementamos el número de Péclet desde cero, los 


tiempos de relajación aumentan, alcanzan un valor máximo y luego disminuyen monóto- 
namente tendiendo a cero cuando el número de Péclet tiende a infinito. 


Asintóticamente se verifica la expresión análoga a la (4.20) del modelo tridimensional: 


li 
2/(t,x)= 3, (1) +) Bio en] (5-0, ] e (4.48) 


u=l 


Si S representa el área de una sección transversal de la cámara y si ubicamos la sección 
de extracción en /—Ó0<x<l, obtenemos unas fórmulas para el modelo global unidi- 


mensional que corresponde a las fórmulas (4.22) del modelo tridimensional: 
t 1 Pe, ( 


” Xx 
A (r; R,) =S- yp efe y a A «Axp4,, (449) 
A partir de las fórmulas (4.49) se puede reescribir la desigualdad (2.11) para el modelo 
global unidimensional en forma análoga a como se hizo en el caso del modelo tridimen- 
sional con la desigualdad (4.23). Al igual que antes, si desde un instante de tiempo en 


adelante se verifica la desigualdad, el coeficiente de enriquecimiento relativo 0, (e ¿R,) 


presenta una fase transitoria por exceso respecto de su valor estacionario 


Mco (R, ) 


19.2 ¿(R.) = —— : en presencia de efectos de difusión cruzados el enriquecimiento 
| 1. KR.) 
1,00 


relativo pasa por un máximo durante el transitorio entre el estado inicial y el estado esta- 
cionario final. 


(5) Solución aproximada del problema de valores y vectores propios me- 
diante la teoría de perturbaciones regulares. Aplicación a las solu- 
ciones de los modelos matemáticos del proceso de enriquecimiento 


Aún en el caso del modelo global unidimensional, los resultados relacionados con la de- 
sigualdad (2.11) resultan difíciles de analizar. 
Tanto los números de moles presentes en la región de extracción en estado estacionario 


n; Ral como la variación en los números de moles respecto del estado estacionario 


Ó. n; (e ,R, ), involucran combinaciones de aportes de los mismos modos normales. Estas 
combinaciones difieren de un soluto a otro solamente en las diferencias que pueden existir 


entre un soluto y otro en las proyecciones e! Y, Un (y = 12. 0) de los vectores 


propios uU A de la matriz de coeficientes de difusión sobre el versor e, característico de 


cada soluto. Estas proyecciones son las i-ésimas componentes de esos vectores propios 


u mE 

Para poder avanzar en el análisis se pueden tomar dos caminos. 

Uno de esos caminos consiste en centrar el análisis en casos concretos, considerando 
soluciones ternarias o cuaternarias ya investigadas experimentalmente y cuyas matrices 
de coeficientes de difusión han sido determinadas para distintos estados de referencia. 
Aplicando métodos numéricos se pueden hallar los valores y vectores propios de las ma- 
trices de difusión (o de dispersión, en el caso del modelo global unidimensional) y luego 
se puede simular el comportamiento de los campos auxiliares a partir de las soluciones 
analíticas halladas en la parte (4), de forma análoga a la empleada en la referencia [2] para 
simular el comportamiento de las concentraciones en el caso particular en el cual los 
efectos de difusión cruzados son despreciables. 

El otro camino consiste en introducir alguna simplificación de carácter más o menos ge- 
neral y continuar la investigación analítica sobre la base de los resultados obtenidos pre- 
viamente en el presente trabajo. 

Este segundo camino es el que seguiremos en lo que resta del artículo. 


En soluciones líquidas y en mezclas de gases diluidos los coeficientes de difusión princi- 
pales suelen presentar valores un orden de magnitud mayor que el de los valores de los 
coeficientes cruzados, aun en soluciones concentradas? [9] [10] [11]. 

Asumiendo la mencionada diferencia en órdenes de magnitud, descompondremos la ma- 
triz de coeficientes de difusión en la suma de una matriz diagonal que presenta en su 


diagonal principal los coeficientes de difusión principales D, li =L 2,...,n) (siendo el 
resto de sus elementos iguales a cero) y una matriz de perturbación formada por los coefi- 
cientes de difusión cruzados: D=D,+e€-D, (5.1) 

Entonces: D, = ls, -D, (5.2) e-D, = ll - 9, ) D, 


Los coeficientes cruzados se escriben D, =E-D, ¡ donde £ es un número positivo de 


(5.3) 


nxn nxn 


orden de magnitud inferior a la unidad y los coeficientes D, ¡¡ Son en, valor absoluto, del 


orden de la unidad o inferior. Como siempre, O, es la delta de Kronecker. 


A partir de una descomposición como la (5.1) es posible aplicar la teoría de perturbacio- 
nes regulares para determinar aproximaciones a los valores y vectores propios de un ope- 
rador [23] [24]. 


3 No obstante, esto puede no ocurrir en el caso de mezclas de gases alejadas de las condiciones de 
idealidad [9] [11]. 


Para ello desarrollamos los vectores propios y los valores propios de la matriz de coefi- 
cientes de difusión hasta el primer orden en € : 


ú, ==, +e-ú, +0,l6) (5.4) D,=D,,+e:D,,+o,(€) (5.5) 
Aquí O ñ (<) es un infinitésimo respecto de £, al igual que la norma del vector O ; (s ): 


Introduciendo (5.4) y (5.5) en D- 7 di Ds -U 7 igualando coeficientes hasta el primer 


orden en € se obtienen dos ecuaciones: D,-u ÓN D 10 -Y dÓ (5.6) 


(5, =D ye 1). Uy e (D,, «I- D)-i,, (5.7) 
En las ecuaciones (5.6) y (5.7) se representa por Í a la matriz unidad. 


Como D, es una matriz diagonal, sus valores propios coinciden con los elementos que 


ocupan la diagonal principal*: Dio =P (u =1 2.0, n) (5.8) 
Sus vectores propios normalizados coinciden con los vectores columna característicos de 
cada componente de la solución: ey coll6,, ¡8 S ” (5.9) 


Consideremos ahora el problema de valores y vectores propios para la matriz D' tras- 
puesta []. 


En general D + D” Entonces D y D' poseen los mismos valores propios D, pero no 
a : E DT > - 
necesariamente los mismos vectores propios: D'-w + D Pd (5.10) 


Los vectores propios de D” también forman un conjunto (7, ,..., *, $ linealmente 


independiente, que junto con el conjunto (a, Md de vectores propios de D se 


pueden normalizar de modo tal que []: w, -Y,= Ono (uo = 1,2....,n) (5.11) 
A su vez la matriz traspuesta se puede descomponer de la siguiente forma: 
D"=D +8sÉ:D (5.12) 


A T » z py T : Ph 
En (5.12) D, =D, mientras que D, es la matriz traspuesta de la matriz D, . 


Introduciendo W,, =W, y +€:*W +0,(€) y D,=D,y+€D, +0,(€) se obtie- 


nl 
ne, al primer orden en € : Woo = Uno (pues D/ = D,) (5.13) 
De (5.13) y de (5.9) se desprende: Wo = (5.14) 
Además, de (5.11) y de los desarrollos en potencias de € : Wo Us Oo (5.15) 


6 Como se dijo antes, supondremos que los coeficientes de difusión principales son distintos dos a 
dos, al igual que los valores propios de la matriz de coeficientes de difusión. Para € lo suficiente- 
mente pequeño, la validez de la segunda afirmación se desprende de la validez de la primera. 


Por su parte Wa verifica: (5, DA 1). Wy = (o, > D,). lo (5.16) 


o . ¿ > T : 
Multiplicando ambos miembros de (5.7) por el vector fila W,,¿ , teniendo en cuenta que 


Wo *D, =D 0 "Wo y la igualdad (5.15), resulta D,¿=Wp*D, y (5.17) 
De (5.14) y (5.3) se desprende que: Dia = e . D, : e, = (u = 1,2,...,n) (5.18) 


— E 7) 7 — — 
Entonces: D,=D, +€:'W,p"D, 4, =D, =D, lu = 12d) (5.19) 
Por tanto, al primer orden en £ resulta que D, = D, 15 cada coeficiente de difusión 
característico de un modo normal resulta igual a uno de los coeficientes de difusión prin- 


cipales de los solutos. 


Al primer orden en € la teoría de perturbaciones regulares [24] suministra las siguientes 


- 


>T > 
2, Wo'D,u 


expresiones para e Y = A z E UL (5.20) 
o=l,0%+u DD, A Do 
Teniendo en cuenta (5.3), (5.14) y (5.19), la fórmula (5.20) se puede reescribir así: 
y o 
7 2 e De, _ a y D, A 
U,=8, €: ya mé, DR e (5.21) 
o=Lo+pu Ds =D y o=lo+u Do —D ip 


Para Wa la teoría de perturbaciones regulares junto con los resultados obtenidos previa- 


: pyT 
mente sobre los valores y los vectores propios de D” : 


=P T 
A das z Uso" D, Woo > ue z Dis > 
EW, = Weo== » EL (5.22) 
o=l,0%pu Doo e* Do o=l,0%u Doo a Din 
Y D 
E . e HO > 
Al primer orden en € : Me E A e DEN (5.23) 


o=l,0%+u Do o Ds 


Para aplicar los resultados obtenidos a las soluciones de los modelos matemáticos tridi- 
mensional y global unidimensional es conveniente introducir un vector columna (vector 
n-dimensional de concentraciones) cuyos elementos son los campos de concentraciones 
de los solutos. Para el modelo tridimensional se tiene: 


cl(t,F)=col(c,(t,F),c,(t,F),...,c, (t,F)) (5.24) 

Para el modelo unidimensional: cl, x) = col(c, (r, x), C), (r, laa (r, x)) (5.25) 
Las concentraciones de los solutos se obtienen proyectando los vectores de concentración: 

c(1,)=87 -¿(0,F) (¡=1,2,....n) (5.26) 

Elt,x)=87 -¿(t,x) (i=1,2,....n) (5.27) 


Los vectores de concentración (5.24) y (5.25) se pueden escribir como combinación lineal 


de los vectores propios de la matriz D utilizando los campos auxiliares como coeficien- 
tes: 


cle,7) =$, (1,7), (528) cle,x)=p,(t,x)- 16, (5.29) 
u=1 u=l 


A partir de (5.28) y (5.29), teniendo en cuenta (5.11) se obtienen los campos auxiliares 
proyectando sobre los vectores propios de D”: 

mw. -clt,r)=9 (t,F) (5.30) mw, -cltx)=0p, (tx) (5631 
Estas últimas fórmulas permiten hallar los valores iniciales de los campos auxiliares a 


partir del valor inicial del vector de concentraciones, cuyas componentes son las concen- 
traciones iniciales de los solutos, iguales en todos los puntos de la cámara. 


Las soluciones tentativas (4.4) y (4.24) se obtienen combinando (5.26) con (5.28) y (5.27) 
con (5.29) respectivamente: 


e (t,7)=8" -¿(1,7)= S 9, (1,7). (2 -5,) (5.32) 
u=1 


a, (t,7)=8) -c(t,x)= y, (t,x)- (er -ú,) (5.33) 


u=1 
Pero de la aproximación (5.21) a U hi al primer orden en £ se desprende que E -Y n es 
igual a 1 si 1=U y si 1% U es de orden de magnitud € , iguala €: A , donde 
ii ua 
e-D, Ss D, " 


Entonces las concentraciones de los solutos para el modelo tridimensional y para el mode- 
lo global unidimensional, al primer orden en £ vienen dadas, respectivamente, por: 


n D. 
e (t,F)= 4, (0,7) + Y 6, (0,7). E (5.34) 


ui, u=1 D, =D. 
clidzoledre: Yoli): Pet — (535) 
pe qe rr Ds =D 


Estos resultados sugieren, para ambos modelos, que el comportamiento temporal y espa- 
cial de los campos C; lc, r) o C; (1, x) es dominado por los campos auxiliares 6, lc, r) O 
O; (e ; x) respectivamente. El aporte del resto de los campos auxiliares aparece como una 


perturbación. 


En (5.35) se explicitaron los coeficientes de difusión efectivos o coeficientes de disper- 
sión, que en el caso de un flujo laminar se relacionan con los coeficientes principales de 
difusión molecular a través de la fórmula (3.19). En (5.34) aparecen directamente los 
coeficientes principales de difusión molecular. 


El número de moles n, 6 ¿R,) del soluto i-ésimo calculados a partir de (5.32) y el núme- 
ro de moles de ese mismo soluto 2, 6 Ey a) calculados a partir de (5.33) vienen dados, 


respectivamente por (por definición [; representa el intervalo O<x<l >): 


Be a e: 2 Zi Dina 
Je (7) avz [6 (1,7). dV+e- Y | 6,(1,7). dv 6530) 


R R, pi, p=1 R, ii uu 


e 


D, 
[c,(t,F). av = S- Jedi) dx+€- Y fo, (t,x): «dx A (5.37) 


l5 WS ii uu 

(Aquí S es el área de una sección transversal de la cámara de enriquecimiento.) 

Las fórmulas obtenidas sugieren nuevamente que el campo auxiliar cuyo índice corres- 
ponde al soluto considerado juega un rol dominante en la cinética del número de moles 
presentes en cada instante en la región de extracción, tanto en el caso del modelo tridi- 
mensional como en el caso del modelo global unidimensional. 


A partir de (5.34) y (5.35) las soluciones (4.14) y (4.46) que describen la evolución tem- 
poral y espacial de las concentraciones de los solutos en la cámara de enriquecimiento se 
pueden expresar en forma aproximada pero más detallada, como sigue: 


(5.38) 
n D. E e t 
N E co rd e Tp. , F 
La DD %. ) 2, Y ) 
E a ES 
G (1, x) PLE) NB 3 E ). 2 sen % )+6,. A 
Ñ (5.39) 


Ahora consideremos los tiempos de relajación Le ó ordenados en forma decreciente: 


E ES E A AA (5.40) 
A partir de (5. 20 se deduce, asintóticamente pea l —> +0 
D, 
mu .. il 
Sn (t,R,)= *-[wal?) )dV+e8- De e A) in (5.41) 
ii uu 
J= 14.0 Vte: » e (5.42) 
TEA 


Reescribamos ió la condición (2.11): si desde un instante de tiempo en adelante 


Sn,(c;R.) _ [6n,(e3R,) 
< 
Mco (R,) m,. (R,) 


presenta una fase transitoria por exceso. Asumiendo la dominancia del primer término 
del miembro de la izquierda de la ecuación (5.41), de la desigualdad (2.11) se desprende 
esta otra, que debe verificarse desde un instante en adelante para que la fase transitoria 
por exceso se produzca: 


entonces la función de enriquecimiento relativo (X,, (e ,R,) 


(5.43) 


Si Tos es menor que T, la desigualdad (5.43) necesariamente se deberá cumplir desde el 


dE 
instante a partir del cual la exponencial decreciente e *'" *'/ se hace inferior al 
miembro de la derecha de la desigualdad. Este miembro no es nulo, puesto que el nume- 
rador es el producto del número de moles del soluto k en la región de extracción, en esta- 
do estacionario por el valor absoluto de la integral ¡AG ) d'V que no se anula. 
Ro 


Una desigualdad análoga se obtiene para el modelo global unidimensional, expresada 


Sn, (e;1,) z Sn, (t;1,) 
a) Moa Ls) 


ahora con referencia la región de extracción /;: 


De esta última desigualdad se desprende esta otra: 


(5.44) 


Pi NS) 


cumplir desde el instante a partir del cual la exponencial decreciente € 
inferior al miembro de la derecha. 


De (4.41) se obtiene la siguiente expresión para el tiempo de relajación dominante del i- 


sde sE l 4. Pe, 

ésimo campo auxiliar: T,=|=|: 3 Z (5.44) 
Pe? +4-77 

Una fórmula igual a ésta se estudió en las referencias [1] y [2], deducida para los tiempos 

de relajación dominantes de los solutos en ausencia de efectos cruzados. Cuando 

Pe, > 277 el soluto cuyo campo auxiliar dominante posee un mayor número de Péclet 

posee un menor tiempo de relajación. 

A su vez, el número de moles del soluto i-ésimo, en estado estacionario y en la región de 

extracción viene dado por (aquí $ es el área de una sección transversal de la cámara): 


n 


Mi, 5110. dire Y 5, A Jo.) dx |b (5.45) 


ui, =1 D, Mu 


Los valores asintóticos de los campos auxiliares vienen dados por la fórmula (4.34), de 
modo que las integrales de esos campos sobre la región de extracción vienen dadas por: 


[0,..2) )-dx=1- to) 20 ox 8] (o =1,2,....n) (5.46) 


Pe 


a 
A integrales (5.46) son funciones crecientes del número de Péclet en el intervalo 
Pe, > 0, que incluye a su vez al intervalo para el cual los tiempos de relajación dismi- 
nuyen cuando ese número crece. 


Ahora bien, por definición (ver (4.28)), Pesa var. Hasta el primer orden en £€ tenemos 


que D, =D,_,,- Ahora, un soluto que posee un coeficiente de dispersión (coeficiente 
principal efectivo de difusión) D, _ menor que todos los demás, posee un valor de 


Pe. que es mayor que el de los demás solutos. 


Si en (5.45) la suma de aportes de primer orden en £ resulte dominada por el término de 
orden cero, cabe esperar que el soluto que posee un mayor número de Péclet se concentre 
más, en estado estacionario, en la región de extracción. 


Para números de Péclet mayores que 277 cabe esperar que este soluto no solamente se 
concentre más en Í, sino que a la vez se aproxime más rápido a su distribución estacio- 


naria, en comparación con los demás solutos. Este resultado justifica la idea de hacer una 
extracción antes de que la distribución de los solutos en la cámara alcance su valor esta- 
cionario. 


(6) Conclusiones 


(a) La condición (2.11) es una condición necesaria y suficiente para que el enrique- 
cimiento relativo ale ,R,) de la solución en el componente k respecto del 


componente Í pase por un valor máximo antes de aproximarse a su valor asintó- 
tico Oz (R, ) A 


(b) En términos de las soluciones de las ecuaciones linealizadas de transporte de ma- 
sa en el marco del modelo tridimensional, la condición (2.11) se reduce a la de- 
sigualdad (4.23). Una desigualdad análoga se puede deducir para las ecuaciones 
de transporte de masa linealizadas del modelo global unidimensional. 

(c) En ambos modelos cada campo de concentraciones de un soluto se expresa como 
combinación lineal del mismo número de campos auxiliares, denominados mo- 
dos normales de convección-difusión, difiriendo un soluto de otro solamente en 
los coeficientes que se utilizan para formar la combinación lineal. 

(d) Cuando los efectos de difusión cruzados son lo bastante débiles, a cada soluto se 
le puede hacer corresponder un campo auxiliar que cumple un rol dominante, 
frente a los demás campos, en cuanto a determinar la cinética de ese soluto en la 
región de extracción. 

(e) En este caso la desigualdad (2.11) se puede reducir a las desigualdades (5.43) pa- 
ra el modelo tridimensional y (5.44) para el modelo global unidimensional. En 
esas desigualdades aparecen el tiempo de relajación dominante del cada campo 
auxiliar correspondiente al soluto considerado. 

(f) Cuando los efectos de difusión cruzados son lo bastante débiles y para números 
de Péclet mayores que 2/7 cabe esperar que el soluto cuyo número de Péclet es 
mayor, no solamente se concentre más en la región de extracción, sino que a la 
vez se aproxime más rápido a su distribución estacionaria, en comparación con 


los demás solutos. 
(g) De la fórmula Pe = EA se desprende que, en principio, es posible aumentar el 


ei 


producto V -Í hasta que los números de Péclet de todos los solutos sean mayores 


que 27 . Como D depende de y (fórmula (3.19)), para aumentar V +! man- 


e,ii 
teniendo constantes los coeficientes de dispersión se debería aumentar la longi- 
tud de la cámara de enriquecimiento. 

(h) Los resultados obtenidos, si bien no justifican un posible interés industrial de la 
idea de hacer una extracción antes de que la distribución de los solutos en la cá- 
mara alcance su valor estacionario, sugieren la posibilidad de un hecho que por 
sí mismo tiene interés científico: la aparición de las fases transitorias por exceso 
previas al estado estacionario. 


En relación con los resultados obtenidos surgen algunas interrogantes: 
¿Si tf =1f,, es el instante en el cual 0, (t ¿R,) alcanza su valor máximo, cuánto vale el 


alt ¡R) 
TES (R, 
ceso de enriquecimiento respecto del proceso original? 


cociente que da una medida del incremento en la eficiencia del nuevo pro- 


¿Es posible seleccionar los parámetros del sistema de forma de obtener una fase transito- 
ria de amplitud y de duración adecuadas para efectuar la correspondiente extracción con 
el máximo de eficacia en el enriquecimiento? 


Para obtener respuestas a estas preguntas se necesita un análisis más detallado, como el 
que se puede realizar efectuando corridas de simulación digital. 

En casos en los cuales los efectos cruzados en la difusión pueden ser despreciados y los 
procesos de transporte se describen mediante un modelo unidimensional, los resultados de 
la simulación digital presentados en las referencias [2] y [3] para números de Péclet ma- 


yores que 2/7 muestran fases transitorias por exceso antes de alcanzar el estado estacio- 
nario. No obstante, la magnitud de los sobrepasos sugiere que la idea de hacer una extrac- 
ción antes de que la distribución de los solutos en la cámara alcance su valor estacionario 
no redundaría en una mejora significativa respecto del proceso de enriquecimiento en 
estado estacionario. 
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